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Abstract 

Analytical solutions for reconnective annihilation in curvilinear geometry are pre- 
sented. These solutions are characterized by current density distributions varying 
both along the radial and azimuthal coordinates. They represent an extension of 
previously proposed models, based on purely radially dependent current densities. 
Possible applications of these solutions to the modeling of solar flares are also dis- 
cussed. 
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1 Introduction 



Magnetic reconnection is a process in which magnetic energy is converted 
into kinetic and thermal energy of a plasma. Its importance lies in the fact 
that magnetic reconnection is believed to be one of the underlying mecha- 
nisms responsible for many phenomena occurring in astrophysical plasmas, 
such as solar flares and geomagnetic substorms [1,2]. After the early works of 
Parker [3], Sweet [4] and Petschek [5] considerable effort has been made to 
improve the theory of magnetic reconnection with the help of exact analyti- 
cal models for reconnective annihilation in Cartesian coordinates [6,7,8]. More 
recently exact solutions for magnetic reconnection in curvilinear coordinates 
have been presented [9,10,11]. Some features of these solutions make them 
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particularly interesting for modeling a large class of solar flares. In this paper 
we present analytical two-dimensional solutions for reconnective annihilation 
in polar coordinates which extend the model derived in [9,10] to the case of 
current density distributions depending on both the radial and the azimuthal 
variable. 



2 Basic equations 



We consider the system of magnetohydrodynamics equations in polar coordi- 
nates for a two-dimensional steady and incompressible plasma with uniform 
density and resistivity discussed in [10]. This system consists of the equation 
of motion 

[^,VV] = [4VM], (1) 

and the Ohm's law 

E+-[rl;,A] = - V V 2 A. (2) 
r 

The flux function A and the stream function ip are related to the radial and 
azimuthal components of the magnetic and velocity field, respectively, in the 
following way: 

In (2) E represents the normalized electric field which is uniform and perpen- 
dicular to the r-9 plane while 77 is the dimensionless resistivity corresponding 
to the inverse magnetic Reynolds number. The Poisson brackets are here de- 
fined as 

if ] = ?L?i - ?i?L (a) 

[J,gl drdO drdO' u 
As shown in [9,10] the ansatz 

A(r,6)=A 1 (r)9 + Mr), (5) 



ij(r,6)=Mr)0 + Mr), (6) 
is compatible with the system of equations (1) and (2). In fact the substitu- 
tion of (5) and (6) into (1) and (2) yields the following set of four ordinary 
differential equations for four unknown functions: 
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- ^Mi + v(Ai' + rA,") = 0, (9) 

E+ l - [^A, + r](A ' + r V)] =0. (10) 

where the dash indicates the derivative with respect to r. In the limit of ideal 
MHD, corresponding to r] — 0, equation (9) implies 

V>1 = CdAi (11) 

with o; being an arbitrary constant. Combining this proportionality relation 
with Eq. (7) yields 

Ax + — ± \ 2 A X = 0, (12) 
r 

where ±A 2 = - (B^dj/dO) I is a parameter related to the variation of the 
current density with respect to the variable 9. For A = the current density 
depends only on r and the solution of (12) is given by 

Ax = ci In r + c 2 . (13) 

For this particular case approximate and exact solutions of the system (7)-(10) 
have been derived (see [10] and [9] respectively). For A ^ the solutions of 
(12) depend on the sign of ±A 2 and they are given by 

A+ = cfJ (\r) + 4Y (\r), (14) 



A- = cr/o(Ar) + c 2 -Xo(Ar). (15) 

where the superscripts of the arbitrary constant C\ and c 2 correspond to the 
sign of ±A 2 . These solutions have an analogous counterpart in the trigonomet- 
ric and hyperbolic sinus solutions derived in [8] for the corresponding problem 
in Cartesian system of coordinates. The ideal solutions for Aq and ipo' can be 
derived from (8) and (10) setting rj = which yields 

,, 1 , a 2 — 1 / b\ t . a Er a b . , 

^o=-^o+ o— \ar+- , A ' = — -— + -r + — , 16 

a a z \ r J a z — 1 A\ a ar 

with a and b arbitrary constants. These solutions are singular in correspon- 
dence to the zeros of Ai, therefore they have an infinite number of singularities 
if the solution Af is considered and only one singularity if ±A 2 ^ 0. It should 
be mentioned that in absence of an electric field the system possesses regular 
solutions describing field-aligned flows with one or an infinite number of null 
and stagnation points depending on the sign of ±A 2 . Singularities in the ideal 
solutions (16) can be resolved by means of a finite resistivity. Here we derive 
solutions for A ^ 0, characterized by current density distributions depending 
on both r and 9, unlike the case A = where the current density is purely 
radially dependent. 
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3 Solutions for A\ and ipi 



For the functions A\ and we assume the following boundary conditions: 

A 1 (1) = B„, ^i(l) = u re , (17) 



A 1 (r c ) = 0, Mr c ) = 0, (18) 

where B rc and t> re represent constant values of the radial component of the 
magnetic and velocity field imposed at the radius r — 1, while r c indicates the 
radius at which Af and ipf have a zero. The ideal solutions (14) and (15), as 
in the case of the logarithmic solution (13), can satisfy all of the four boundary 
conditions (17) and (18) taking the form 

A+ = B J (\r c )Y (\r)-Y (\r c )J (\r) 

rc Y (X)J (\r c )-Y (Xr c )J (X) { } 

and 

a- _ t> Io(\r c )K (\r)-K (\r c )I (\r) 
Al ^ c K (X)I (Xr c ) - K (Xr c )I (X) ' ^ 



respectively, together with 



^ = (21) 

-Ore 



A fundamental difference between the cases A = and A 7^ is that for the 
former the ideal solutions for A\ and ip\ provided by (11) and (13) are also 
exact solutions of the complete resistive equations (7) and (10). For A ^ 0, 
however, this is not the case. In fact, comparing the Taylor expansion about 
r c of the exact solutions of (7) and (9) with the one of the ideal solutions it 
transpires that a mismatch occurs between Af"(r c ) and the third derivative 
of the exact solution for A 1 at r c . Due to the high complexity of the nonlinear 
equations (7) and (9) we limit ourselves to solve them in an approximated way. 
In particular we seek solutions satisfying the imposed boundary conditions and 
whose Taylor expansions about r c are compatible with the ones of the exact 
solutions at least up to the third derivative, in order to resolve the lowest 
order mismatch between ideal and exact solutions. Since the typical value of 
the dimensionless resistivity r] for many astrophysical and laboratory plasmas 
is between 1CT 10 and 1CT 14 we can consider i] as a small parameter and use the 
method of matched asymptotic expansions. We separate the problem in two 
domains, an inner region consisting of a narrow layer containing r c , where the 
ideal solution breaks down, and an outer region, complementary to the inner 
region, where (19) and (20) are approximate solutions. Asymptotic expansions 
for Ai and ipi are found in the inner region and then matched with the outer 
zeroth order expansion (19), (20) and (21). In order to derive inner solutions of 
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the nonlinear equations (7) and (9) the dependent and independent variables 
are rescaled in the following way: 



r-r c . 



(22) 



The equations obtained after the rescaling are 
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It should be noted that in the limit r c — > oo the above equations reduce to the 
rescaled equations for the problem in Cartesian coordinates considered in [8]. 
An asymptotic expansion in powers of r] is chosen as a form of solutions for (23) 
and (24). Also in curvilinear geometry the rescaled equations depend explicitly 
on i] and in particular the square root of i] is present in them. Therefore it is 
natural to choose inner solutions of the following form: 



M = A 10 (s) + y/rjA u (s) + r]A 12 (s) + • 
ipi = ipio(s) + y/rjil>ii(s) + rji) l2 {s) + • • 



(25) 



where, unlike to the case of Cartesian geometry, the square root of rj appears 
in the expansions. The resulting inner solutions are then written in terms 
of the original variable, expanded as rj — > for fixed r and matched with 
the expansion about r c of the outer solutions (19) and (20). Eventually one 
obtains 
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^re^c 

is a positive quantity. The inner solutions (26) possess the required properties, 
since they satisfy the boundary conditions (18) and Ai ± '"(r c ) is consistent with 
the value of the third derivative of the exact solution at r c . By combining 
(26) with (19), (20) and (21) composite solutions can be found, which are 
approximately valid over the whole domain of interest. 



4 Solutions for A ' and ip ' 



As boundary condition for A ' we impose 

A '(r c ) = 0. (29) 

This constraint, together with (18), fixes the position of a magnetic null at 
the point (r = r c ,9 = 0). On the other hand, as mentioned in Sec. 3 the 
ideal solutions (16) become singular at r c which is in striking contrast with 
the condition (29). Therefore the ideal solution for Aq and -?/V also breaks 
down at r = r c and a boundary layer analysis is required in order to find 
approximate resistive inner solutions satisfying (29). The inner equations for 
Aq and ipo are derived by rescaling the independent variable in Eqs. (8) and 
(10) in the following way: 



r - r c 




Then an asymptotic expansion in powers of r] is sought in the form 

A ' = < + ^jAl', Vo' = ^o°' + Vv^o- (30) 

In the equations for the inner solutions A\ and ipi are replaced by the leading 
terms of the Taylor expansions of (26) about s = 0. The resulting equations 
are formally equivalent to the corresponding equations solved in [10]. The 
composite solutions obtained after the matching with (16) are given by 



6 



A ' = E 



daw(v/c ± s 



2 1 

s + 



3 r c rev^ Vk±r c 2r c fc ± y/rj\ k± 



1 , ±s . fc± 



+ — s - 



+ 



+ 



r c fc ± V3 3 6^ k±^,s A;±r c 2 lnf^s_ 



- \ar c + 



b \ B Te k ± s 



+ 



a(sJ2rj + r c ) + 



= Aq + 

Vre 



ar c H . 



5 Discussion of the solutions 



The solutions derived in Sec. 3 and 4 describe stationary magnetic and veloc- 
ity fields whose topologies are nontrivial due to the presence of a magnetic 
null point at (r = r c , 6 = 0) and of a stagnation point at (r = r c , 6 = 
[(— v rc 2 — B IC 2 )/v TC 3 a ± ](ar c + b/r c ). Each of these singular points is located 
at the intersection of two separatrix lines, one of which, corresponding to the 
arc r = r c , is in common to both points. Along this shared separatrice a 
curved sheet of enhanced current density is present, allowing for reconnection 
of magnetic field lines that are advected by the flow from the outer region. 
The rate at which reconnection occurs is given by the parameter E that does 
not depend on rj. Thus these solutions belong to the class of fast reconnection 
models. As in the case of many other magnetohydrodynamics models of re- 
connective annihilation [7,8,9,12] a limit to the reconnection rate is imposed 
by the requirement of incompressibility. Indeed, consider the expression for 
the plasma pressure given by 

v r 2 B e 2 [ B e 2 - 

with constant po- The second term on the right hand side of (31) is negative 
and diverges as r — > oo. Therefore, in order for p to be a positive quantity, 
the domain of the solutions must be limited to a finite region of space. 
The current density is perpendicular to the r-6 plane and its expression is 
given by 

i = ~ (t + A ") e ~ V " A ""- < 32 > 

from where it is possible to notice the linear dependence of j on 9. With the 
help of (12) it is also possible to see that this dependence is absent for A = 0. 

The variation of the current density with 9 can clearly be seen in Fig. 1, which 
shows that this variation is considerably large in the outer region, while in the 
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Fig. 1. Magnetic field lines (thick) and streamlines (thin) for A = 8.7, E = 0.5, 
r) = 10" 2 , v ie = 0.8, B vc = 0.9, r c = 0.2, a = and b = 1. The plots of the field 
lines are superimposed to the distribution of the current density in gray half-tones. 

inner region the current density is practically only radially dependent. 
Finally, from Fig. 1 one can also notice that for the branch of solutions match- 
ing A\ + and the resulting magnetic configuration can be considered as 
generated by an odd number of sources with alternating polarities lying on 
the plane y = 0. The number of sources depends on the value of A and on how 
large the domain of validity of the solution is. The latter cannot be chosen 
arbitrarily large because it would then include additional zeros of (14) where 
the solutions become singular. However this is not an additional restriction. 
Indeed, as already mentioned above, the exact magnetohydrodynamic solu- 
tions for incompressible plasma are usually valid only in a finite size region 
anyway because of the lower limitation on the pressure. On the other hand, 
the magnetic and flow configurations corresponding to ±A 2 < resemble very 
much the ones for A = solutions, which are regular everywhere and therefore 
they are not subject to restrictions due to the presence of singularities. The 
current density for this branch of solutions varies also linearly with 9 but in 
the opposite direction to the case ±A 2 > 0. 



6 Conclusions 

Solutions of the steady two dimensional MHD equations for an incompressible 
plasma in polar coordinates have been derived using a matched asymptotic 
expansion technique where the inverse magnetic Reynolds number has been 
assumed to be very small. Two branches of solutions have been obtained de- 
pending on the value of one parameter related to the variation of the cur- 
rent density with the azimuthal variable. The magnetic field configuration 
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described by the solutions has a null point at a certain radius. Along this 
radius a curved current layer is formed, which causes the reconnection of field 
lines. In the magnetic configuration three alternating polarities lying on the 
plane y = produce the magnetic field with the null point in the volume 
above. This feature has interesting applications to solar physics because the 
presence of three alternating polarities on the photospheric plane is a feature 
common to a large class of solar flares [13]. For one branch of solutions, in 
general, on the plane y = an odd number of sources is present, three of 
which interacting with each other. This number depends on the value of the 
parameter A and on the size of the domain of validity of the solutions. The 
reconnection rate does not depend on the resistivity therefore these solutions 
belong to the class of models providing fast reconnection regimes. 
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